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Abstract. The kinematics of distant OB stars perpen- 
dicular to the Galactic plane, inferred from proper mo- 
tions in the Hipparcos catalogue, are analysed and com- 
pared to the kinematic signature that would be induced 
by a long-lived Galactic warp. Previously we reported 
that the kinematics of the OB stars toward the anticen- 
ter were inconsistent with the assumption of a long-lived 
warp ( fgmartetal. 1998 ), showing negative systematic mo- 
tions as opposed to the expected positive motions. Using 
a larger sample of OB stars, improved distances, and a 
more accurate model for the uncertainties, we confirm our 
previous result for a larger range of galactocentric radii. 

However, we note that the new model for errors in 
the photometric distances reveal an important bias that 
causes the observed systematic vertical motions to be 
smaller than their true values. Using synthetic catalogues 
we investigate the effect of this bias on the systematic 
vertical motions in conjunction with the possibility of a 
smaller warp amplitude and a warp precession rate, both 
of which can also lead to smaller systematic motions. 
Taken together these three effects can work to produce 
negative observed systematic vertical motions, similar to 
those detected in the data, though only with both exces- 
sively high precession rates (— 25 kms -1 kpc -1 ) and very 
large photometric errors (1 magnitude). 
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1. Introduction 

For some time now the existence of a warp in 
our Galaxy has been well known. First seen in 
neutral hydrogen emission, it's presence in sev- 
eral other components have since been reported, 
i.e. IRAS point sources (pjorgovski fc Sosin 1989^ , 
CO gas (Woutcrloot et al. 199C), and dust 



(Freudenreichet al. 1994). Studies of the Galactic warp 
have been mostly limited to it's spatial structure (see 
Binncy 1992 for review) because the kinematic signature 
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of a warp primarily manifests itself in the component 
tangential to the linc-of-sight (LOS) of an observer 
near the Galactic plane, the component that cannot be 
directly measured by radio observations. For stars, on 
the other hand, this component is in principle accessible 
via proper motions, but in practice this approach has 
been limited by the small magnitude of the expected 
systematic motion (~ 1-2 masyr -1 ), which is smaller 
than the zonal systematic errors in traditional catalogues. 
The advent of astrometry from space, inaugurated with 
the Hipparcos satellite, signals a new era in the study of 
Galactic kinematics; the high precision (< 1 masyr -1 ) 
and accuracy (< .1 masyr -1 ) of the proper motions 
combined with and inertial frame means that it is now 
possible to detect the kinematic signature expected from 
a Galactic warp. Once the spatial structure of the warp 
is determined, its kinematic signature will yield direct 
information on the precession rate of the warp. 

Earlier we reported results from a subsample of the 
Hipparcos catalogue ( Bmart et al. 1998 ). In that study we 
looked exclusively at distant OB stars toward the anti- 
center (between 70 and 290 degrees in galactic longitude) , 
expecting such stars to trace the motions of the gaseous 
component from which they were recently born. Unex- 
pectedly it was found that their kinematics did not follow 
the predicted signature of a long-lived warp, either pro- 
cessing or not. This led to the conclusion that either there 
are additional systematic vertical motions in the Galaxy's 
gaseous component, or that the warp is not a long-lived 
structure. 

In this contribution we use a larger sample of OB stars 
from the Hipparcos catalogue to reinvestigate their sys- 
tematic vertical motions. As before, our interpretation of 
the data is based upon a comparison between the ob- 
served kinematics of the selected sample of OB stars, 
and the "observe" kinematics of a sample created using 
a synthetic catalogue generator. The synthetic catalogues 
are generated based on an assumed statistical model of 
the (warped) stellar distribution, kinematics and observa- 
tional errors. By making the identical selections and cuts 
on the simulated catalogues as are made on the actual 
data, any biases present or introduced in the sample of 
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OB stars due to errors or selection criterion are likewise 
reproduced in the synthetic samples as well. In this way 
the modeled kinematics are translated into observed quan- 
tities which can be directly compared with those of the 
selected sample of stars. In essence, rather than attempt- 
ing to remove or correct biases in the observed sample, 
we introduce the same biases into the model in order to 
determine what is actually detectable. 

In the following section the new sample of OB stars is 
described, while in Sect. 3 the model for the errors used 
for the synthetic catalogues is presented. Sect. 4 details 
an improved estimate of the distances and in Sect. 5, us- 
ing a simple model for the Galactic warp, we find spatial 
parameters consistent the distribution of the Hipparcos 
OB stars. In Sect. 6 we discuss the estimation of the sys- 
tematic vertical motions, derive the kinematic signature 
resulting from a long-lived precessing warp, and compare 
the observed and expected kinematic signatures. Sect. 7 
describes the bias introduced in the estimation of the ver- 
tical motions by errors in the distance and how these, 
in conjunction with a smaller amplitude and precession, 
can produce smaller systematic motions than otherwise 
expected. We present our conclusions in Sect. 8. 

2. The Sample 

The sample of OB stars used here to study the Galactic 
warp is taken from the 10544 OB stars in the Hipparcos 
catalogue; after excluding 3 stars whose spectral typing 
is no more specific than "O" or "B" , we find 4538 stars 
with a parallax 7r < 2 mas. Stars with measured parallaxes 
7r < are retained, as their exclusion would unfavorably 
bias the sample with respect to distance. 

In the Tycho catalogue there are 857 stars brighter 
than V magnitude 7.5 which do not appear in the Hip- 
parcos catalogue, which has a total of 24384 stars brighter 
than this magnitude. This suggests that the Hipparcos 
catalogue is approximately 97 percent complete to this 
magnitude limit. From our sample alone, based on the 
distribution of apparent magnitudes brighter than magni- 
tude 7, we estimate that the sub-sample of 3840 OB stars 
with m < 7.5 is 98 percent complete ( [Smart ctal. 1999| ). 
Of these brighter stars, 929 have a parallax tt < 2 mas. 
Fig. [j] shows the distribution of this subsample of distant 
OB stars on the sky. 



In our earlier work ( Smart etal. 1998) purely photo- 
metric distances were used. In Sect. 4 we will improve 
upon these distance estimates, but nevertheless complete 
spectral classification is needed to calculate a photomet- 
ric parallax. For the 7279 OB stars with complete spectral 
classification supplied by Hipparcos, the absolute magni- 
tudes and colors from Schmidt-Kaler (1982) are used to 
calculate the photometric parallax. For 141 stars without 
luminosity classification H/3 linewidths were used to de- 
termine luminosity class. The 3121 remaining stars with- 
out luminosity classification have photometric distances 
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Fig. 1. Distribution of the subsample of 929 distant 
(7r < 2 mas) Hipparcos OB stars on the sky with apparent 
magnitudes brighter than 7.5 . 



derived from their reddening and estimated absorptions. 
For further details of these procedures the reader should 
consult Smart et al. (1997). 

In order to understand the properties of the photomet- 
ric error, the sample of all OB stars brighter than 7.5 will 
be considered in the following section. These are 3840 in 
number, with 789 not having an initial luminosity classifi- 
cation (nonLC stars). As we expect that our estimation of 
the luminosity classes for these stars will be more suscep- 
tible to error than those with full spectral classification, 
it is important to be aware of biases in their distribution. 
Indeed, because of the us e of the Michigan Spectral Sur- 
vey in the southern sky (Houk fc Smith-Moore 1988 and 



references therein), 745 of the bright nonLC stars (94 %) 
are in the north (5 > —12. deg). This zonal bias results 
in a photometric error that is systematically larger in the 
northern part of the sky than the south, as the following 
section describes. 
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Fig. 2. The distribution of nonLC stars in the north (S > 
— 12. deg) with respect to spectral type. 
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Fig. 3. The distribution of the northern nonLC stars in 
apparent magnitude for the spectral types B3, B5, B8 and 
B9. 



To account for the higher photometric error introduced 
by the nonLC stars, it is necessary to describe their dis- 
tribution not only on the sky, but also with respect to 
other relevant variables. Fig. shows the relative fraction 
of northern nonLC stars with respect to spectral type and 
Fig. |^ the relative fraction of nonLC stars with respect to 
apparent magnitude for B9, B8, B5 and B3 spectral types. 
The fact that the nonLC stars are primarily restricted to 
these types, and not hardly seen in types B7, B6, and B4, 
indicates that the nonLC stars are in general typed at 
a lower resolution than those stars with luminosity class. 
For stars of spectral types earlier than B2 the relative frac- 
tion of nonLC stars does not vary strongly with apparent 
magnitude. 

While in Sect. 4 only a magnitude cut is applied, in 
Sects. 5 and 6 two additional cuts are employed; stars with 
a parallax of it > 2 mas and a vertical velocity greater 
than 50 kms -1 are removed. The parallax cut is im- 
posed to avoid being biased by local structure, such as 
Gould's Belt, which dominates the distribution of nearby 
OB stars. The velocity cut is astrophysically motivated 
by the presence of runaway stars, which through close en- 
counters have reached escape velocities, and thus do not 
possess kinematics principally determined by the Galactic 
potential. The typical threshold detection value for the rel- 
ative space velocity of runaway stars varies from 40 to 65 
kms -1 (Blaauw 1993, Torra et al. 1997), thus our value 
of 50 kms -1 applied to one component is conservative, 
especially as it is approximately seven times greater than 
the velocity dispersion of this population. The velocity cut 
is applied to galactocentric radial bins, excluding those 
stars which deviate from each bin's median velocity by 
more than 50kms _1 , reducing the sample of bright, dis- 
tant OB stars to 895 in number (less 34), 142 of which 



are nonLC stars. Meanwhile 255 are cut from all the stars 
with 7r < 2, leaving 4283 stars, 1291 of which are nonLC 
stars. The calculation of the individual and mean vertical 
velocity is discussed in Sect. 6. 

3. Simulated catalogues and improved error 
model 

In order to interpret the observed kinematics of OB stars 
from the Hipparcos catalogue, a program was developed 
to generate synthetic catalogues based on a statistical 
model of the spatial and velocity distributions of the OB 
stars (Drimmel et al. 1997, 1999). The model distribution 
is non-axisymmetric, being described by spiral arm seg- 
ments, most important of which is the local Orion arm, 
which has a pitch angle of 8 degrees and a spur projecting 
toward the outer Perseus spiral arm. The synthetic cata- 
logue generator produces catalogues with distributions in 
galactic longitude and latitude (I, 6), apparent magnitude, 
and proper motion, which closely resemble those of the 
distant (ir < 2 mas) OB stars in the Hipparcos catalogue 
(see also (Smart ct al. 1997 ) . 

The synthetic catalogue generator uses a Galactic dust 
distribution model based on the 240/j, COBE data (Drim- 
mel & Spergel, in preparation) to calculate absorbtions. 
This model is a refined version of the earlier work of 
Spergel et al. (1997), adding nonaxisymmetric structure 
to the dust distribution model, the most important of 
which is dust associated with the local Orion arm. The 
added structure in the dust distribution results in addi- 
tional reddening in the directions tangent to the local arm, 
thus improving the modeled longitude distributions (see 
Drimmel ct al. 1999 as compared to Drimmel et al. 1997). 
This same model is used to estimate reddenings and ab- 
sorptions to the nonLC stars. 

Observational errors are modeled for generating ob- 
served quantities. The simulated catalogues in Smart et al. 
(1998) used parallax errors, Ot, equal to 1 mas, and proper 
motion errors, 07, , were described by simple functions 
of galactic latitude, as seen in the Hipparcos catalogue 



(Mignard 1997). These errors are actually better described 
as varying with ecliptic latitude and apparent magnitude. 
We have thus improved our model for <7 t and by inter- 
polating from Tables 3.2.4-6 given in the first volume of 



the Hipparcos catalogue (ESA 1997). We have also added 
a 0.01 magnitude error to the apparent magnitudes. For 
our purposes here the Hipparcos errors in position (I, b) 
are negligible. 

However, the improvement of the model of photometric 
distance error makes the greatest difference in the simu- 
lated observable kinematics of the synthetic catalogues. 
Previously our error model involved adding a relative 
gaussian error of 25 percent of the distance. However, if 
the errors are normally distributed in absolute magnitude, 
it is the errors in the distance modulus that are normally 
distributed. We now model the error in photometric dis- 
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Fig. 4. Photometric versus trigonometric parallax distributions of the bright (m < 7.5) Hipparcos OB stars for the 
north (S > —12. deg) and south. The dashed line has a slope of 1, given as reference, while the solid line is the result 
of a robust linear fit that minimizes the absolute deviations (no points excluded). 




Fig. 5. The photometric-trigonometric parallax distributions of the bright Hipparcos OB stars in the north, decom- 
posed into stars with initial luminosity classes (left plot) and without initial luminosity classes (right plot). Dashed 
and solid lines are as in the previous figure. 



tance by assuming a normal distribution of errors in the 
absolute magnitude. This error is not only the result of the 
intrinsic scatter of absolute magnitudes, termed "cosmic 
error" , but also due to errors introduced by uncertainties 
in the determination of the absorption, intrinsic scatter in 
the colors, and misclassification errors. Not all these errors 
are guaranteed to be gaussian. 

The photometric error of our observed sample of stars 
is evaluated using the photometric parallax np verses 
trigonometric parallax 7r t distribution for the sample of 
3840 OB stars brighter than 7.5. Fig. |I] shows the 7rp-7r t 
distributions for the north and the south. The degraded 
photometric parallaxes of the north, due to the presence 



nonLC stars, are responsible for the larger scatter as com- 
pared to the south. To illustrate this effect Fig. || shows 
the 7Tp-7Tt diagrams for the north decomposed into the 
stars with and without luminosity class. It can be seen 
that the northern stars which possess full initial classifi- 
cation (Fig. 5 left) have approximately the same scatter 
as the southern stars (Fig. 4 right), and we will assume 
the same error model can be used for both. One difference 
between the northern and southern stars with full initial 
classification is a systematic bias of 7Tp in the south for 
stars with n t centered at 7.5 mas. This bias is due to the 
presence of stars in the direction of p Ophiuchus whose 
absorptions are systematically overestimated. 
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Fig. 6. Modeled photometric-trigonometric parallax distribution of bright OB stars in the north, for LC and nonLC 
stars (left and right plots respectively), as produced by the synthetic catalogue generator. Dashed and solid lines are 
as in the previous figures. 



Assuming that the observational errors in the Hippar- 
cos trigonometric parallax are well described, we can esti- 
mate the photometric error using the synthetic catalogue 
generator, comparing observed and generated 7rp-7r t dis- 
tributions. For the northern stars with luminosity classi- 
fication an error model which approximately reproduces 
the observed distribution of trigonometric and photomet- 
ric parallaxes is 

a M = .5 + max(0, (to - 6.5)/6) - M/12, (1) 

where the apparent magnitude term can be interpreted 
as being the result of a higher occurrence of classifica- 
tion error at fainter magnitudes, and the absolute mag- 
nitude term can be interpreted as an increasing cosmic 
error and/or uncertainty for the more luminous stars. To 
characterize the distribution a robust linear fit is made 
to the observed np-nt distribution in the north (Fig. |J), 
giving a slope of 0.77 and a mean absolute deviation in 7r t 
of 1.28, after 17 stars are excluded whose absolute devi- 
ations in 7Tt are greater than 6 mean absolute deviations 
(presumably misclassified stars). Meanwhile the modeled 
distribution (Fig. ||) produces a slope of 0.80, with 1.06 
being the mean absolute deviation. For comparison, Fig. 

shows the 7rp-7r t distribution for a model generated with 
a purely gaussian magnitude error of <tm = 1 magnitude, 
showing that such a high error does not reproduce well 
the observed 7rp-7Tt distribution. 

For the nonLC stars the signatures due to misclassifi- 
cation are prominent, its nongaussian nature leaving clear 
artifacts in the distribution, especially at n t > 2 mas and 
7rp < 3 mas (Fig. Bright). These are intrinsically dim stars 
incorrectly assigned a high luminosity. In contrast, intrin- 
sically bright stars misclassified as low luminosity stars do 
not produce clear signatures, but do cause stars of low 7Tt 
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Fig. 7. Modeled photometric-trigonometric distribution 
of bright OB stars in the north assuming um — 1 magni- 
tude. Dashed and solid lines are as in previous figures. 



to be systematically biased toward higher 7Tp. We model 
the misclassification errors as offsets to the true absolute 
magnitudes, with a binned gaussian probability of mis- 
classification. 

However, before applying such errors the distribution 
of nonLC stars must be modeled. The synthetic catalogue 
generator does not at this time generate colors, so spec- 
tral types were simply correlated with the main sequence 
absolute magnitudes to infer a model of the distribution 
of nonLC stars in the north. Fig. || shows the model of the 
relative fraction of nonLC stars in the north. Once a star 
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1 .o Table 1. Modeled offsets in absolute magnitude due to 

misclassifications, used according to value of the binning 




Fig. 8. Model of the relative fraction of nonLC stars in the 
northern part of the sky, as used in the synthetic catalogue 
generator. 

is randomly determined to be nonLC, a binning error is 
calculated using 



p = int(1.5iV + 31ogd + 1-4), 



(2) 



where TV is a gaussian deviate of unit variance and do the 
true distance, (31ogdo + 1.4) being a bias reflecting that 
nearby stars are more likely to be given an over-luminous 
classification, and distant stars a under-luminous classifi- 
cation. The resulting index p is then used to indicate the 
appropriate offset, as given in Table ^. In addition to this 
offset, gaussian noise of <jm — 0.3 magnitude is added to 
the absolute magnitudes. 

Much of the added complexity of misclassification in 
the nonLC stars will not be needed in what follows, as 
for stars with 7r < 2 mas only the positive offsets will con- 
tribute significantly to the error. 

4. Improved distances 

In order to improve upon the distances, a distance esti- 
mate is used that corresponds to the weighted mean par- 
allax of the trigonometric and photometric parallaxes: 

1 TTt/^t + TTp/d 

-, = TTw = 

a 



IK 



l/a P 



(3) 



where ap — 7rpOA//2.17 (Smith 1985), and at is interpo- 
lated from Tables 3.2.4-6 in volume one of the Hipparcos 
catalogue. For this distance estimate stars with initial full 
spectral classification are assumed to have a <jm as spec- 
ified by Eq. (|l]), while nonLC stars are assigned <jm = 3 
magnitudes. This prescription is used for both synthetic 
and observed catalogues. 

To test whether this new distance estimator is more or 
less biased than the photometric distance alone, the syn- 
thetic catalogue generated in the previous section for the 



northern sky is used to compare the estimated distance 
with the true distance. Figs. || and |l(] show the result of 
this Monte Carlo test; for stars with an initial luminos- 
ity class, the variance and skew of the deviations from 
a robust linear fit to the photometric distances, are 2.02 
and 1.04, while with the new distance estimate they are 
0.41 and -0.07 respectively, showing that the uncertainties 
are effectively reduced. The slopes and intercepts of both 
distance estimates show no significant systematic biases. 
For the nonLC stars the new distance estimator reduces 
the scatter produced by misclassifications, though the dis- 
tances can still have significant error. 

Using the typical errors <jm — -7 and at = 1.2 mas, the 
error of 7r w for a star at 7r t = 7rp = 1 mas is about .3 mas, 
for which the distance estimate would be lt^kpc. We 
also point out that for stars in our irt < 2 mas sample, 
the photometric distance will dominate in Eq. (|Io|), as it 
is only at a distance of d = a m / '(2.17 at) that ap = a t . 

5. Warped spatial distribution 

As the Sun is near the line of nodes of the Galactic warp 
the vertical deviation from the Galactic plane is small, but 
a local slope with respect to the Galactic plane is produced 
in the stellar distribution. Using the subsample of the 894 
bright (m < 7.5 magnitudes), distant (nt < 2 mas) OB 
stars, the slopes of the stellar distribution are found in 
the heliocentric coordinates z and y, y being taken in the 
direction of rotation, for bins in galactocentric radius r, 
(r© = 8 kpc) . The slopes were found using a standard 
robust linear fitting routine which minimizes the absolute 
deviations in y. In Table |^ we express these slopes as a tilt 
angle 9, corresponding to the tilt of a galactocentric ring 
whose axis of tilt goes through the Sun, consistent with 
the assumption that the line of nodes of the Galactic warp 
goes through the Sun. 

It is well known that OB stars are preferentially found 
in associations, resulting in a "dumpiness" of the distri- 
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Fig. 9. True versus weighted mean distance (Eq. (|3|)) of stars in a synthetic catalogue for stars with initial luminosity 
class (left plot) and nonLC stars (right plot). Dashed and solid lines are as in previous plots. 
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Fig. 10. True versus photometric distance of stars in a synthetic catalogue for stars with initial luminosity class 
(left plot), and nonLC stars (right plot) where misclassification artifacts are evident. Dashed and solid lines are as in 
previous plots. 



bution on a larger scale than would be found in a kinemat- 
ically relaxed population. Nevertheless, we do not expect 
the presence of associations to affect our results above, as 
we are describing the overall distribution of OB stars on 
a scale significantly larger than the size of a typical asso- 
ciation. To test this presumption stars which were found 
to be within the spatial limits of a given association were 
identified. For this purpose the positions (I, b, d) and limits 
(±AZ, ±A6) for 38 associations were taken from the sum- 
mary list in Lang (1992), while the line-of-sight extent 
of each association (Ad = dmax(AZ, A6)) was doubled to 
take into account possible errors in the distance. Only four 
associations have ten or more stars from our complete sub- 
sample, these being the Cygnus 7 (10 stars), Carina (13 



stars), Cepheus 2 (19 stars), and Orion (17 stars) asso- 
ciations, where a total of 102 stars are found inside all 
associations. Removing each of these sets of stars in turn, 
the tilts were redetermined. It was found that only one 
of the associations (Cepheus 2) influenced the measured 
slopes by more than 0.5 deg. 

A warp is effected in the modeled distribution by de- 
scribing the vertical distribution as a function of z' = 
z — Z w (r,(j}), where the function Z w (r,(f>) describes the 
spatial form of the warp in galactocentric cylindrical co- 
ordinates, <p being taken in the direction of rotation. From 
radio observations the warp is seen to be symmetric out 
to about r — 16 kpc, and can generally be described by 

Z w — h(r) sm(4> - <j) w + uj p t), (4) 
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Table 2. Angular tilts, 6, of galactocentric rings, in units of degrees, implied by the distribution of bright (m < 7.5), 
distant Hipparcos OB stars in heliocentric z and y coordinates. (9 — arctanfr, where z = a + by) For the warped 
and non-warped cases (columns 4, 5, and 6) the tilts and uncertainties are the means and standard deviations of 30 
simulated catalogues. The warped cases (columns 5 and 6) were generated with r w — 6.5 kpc and = 15 and 30kpc 
respectively (see Eq. (|4|)). Column 2 gives the number of stars in each bin for the data; the 92 stars not appearing in 
the table are outside the range of the bins. 



Bin (kpc) 


N 


Data 


No warp 


Warp 


1/2 Warp 


6.5< r <7.5 


140 


-0.39 


-0.24 ± 0.94 


0.34 ± 0.91 


0.29 ± 0.73 


7.5< r <8.5 


500 


1.67 


-0.34 ± 0.44 


1.08 ± 0.49 


0.25 ± 0.44 


8.5< r <9.5 


161 


2.05 


0.07 ± 0.93 


2.54 ± 0.86 


1.45 ± 0.80 



4> w being the phase of the warp, and lo p the precession rate 
of the warp in the direction opposite to Galactic rotation, 
and with 



h(r) = 



(r 
0. 



r > r v 
r <r v 



(5) 



being a height function parameterized by and the ra- 
dius at which the warp starts, r w . The radio observations 
show that the Sun lies close to the line of nodes, that is, 
\4>w\ < 10 degrees. As th e systematic velocity is n ot sensi- 
tive to <f> w in this range ( Smart fc Lattanzi 1996 ) we have 
assumed (f> w = in this study. The assumption of a long- 
lived warp is effectively implemented by taking lo p to be 
a constant with respect to r and time t. 

Table || shows the mean tilts obtained from thirty syn- 
thetic catalogues with and without a warp, one warped 
distribution being generated with the parameters — 
15 kpc and r w = 6.5 kpc, while a second model with half 
this amplitude (r^ = 30 kpc) is also shown for comparison. 
The first set of parameters produce tilts that are consistent 
with the observed tilts, and these are adopted as our esti- 
mated warp parameters in what follows. As noted in our 
earlier work, the observed tilts indicate a Galactic warp 
starting within the Solar Circle. 

The use of a robust linear fitting routine does not pro- 
vide formal errors on individual fits, hence no error is given 
on the tilts of the data. In the case of the synthetic cat- 
alogues we have the freedom to directly evaluate the un- 
certainty by generating multiple catalogues, as we have 
done. Inasmuch as the synthetic catalogues are accurate 
statistical representations of the observed data, the uncer- 
tainties of their tilts will be representative of the uncer- 
tainties in the tilts of the observed catalogue. However, the 
synthetic catalogues do not include the dumpiness due to 
the presence of OB associations, which will increase the 
uncertainties. The 0.5 deg influence of the Cepheus 2 as- 
sociation on the tilts can be taken as an estimate of the 
probable additional random error that should be added 
to the uncertainties quoted in Table 0. In any case, the 



observed tilts are significant and inconsistent with the no- 
warp hypothesis. 

Having estimated the parameters that describe the 
spatial form of the warp in the OB stars, we are now pre- 
pared to consider the kinematic signature accompanying 
such a warp. 

6. Kinematic signature of the warp 

The presence of a long-lived warp in the stellar disk will 
result in systematic motions perpendicular to the Galac- 
tic plane; for an observer in the plane of the Galaxy 
such motions will be primarily observable in the galac- 
tic latitude component of the proper motion. Taking x = 
(cos I cos b, sin / cos &, sin b) as the unit pointing vector to a 
star, the unit vector tangential to the line-of-sight (LOS) 
in the direction of increasing galactic latitude is 



x x (k x x) k — x sin b 



|x x (kxx)| 



cos b 



(6) 



where k = (0,0,1), the unit vector perpendicular to the 
Galactic plane. The component of a star's relative motion 
in the direction of b is 



(v - v Q ) • b = 4.74dfi b , 



(7) 



where v Q is the solar motion, v the star's velocity with 
respect to the local standard of rest (LSR), d and //& being 
the distance in kiloparsecs and the galactic latitude proper 
motion in masyr -1 respectively. Taking v = (U, V, W) Q 
and v = (U, V, W), defined in the same right hand coor- 
dinate system as x, one can write from Eqs. (^) and (^): 



W 



A.74d/J, b 
cos b 



+ W Q + (S - So) tan 6, 



(8) 



where S is the component of the star's velocity parallel 
to the Galactic plane and in the plane that contains the 
LOS and is at right angles to the Galactic plane, that is 
S = Ucosl + Vs'ml, and similarly for Sq. Since distant 
stars are being considered, we do not wish to assume that 
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the stars are in the solar neighborhood (Oort's approxi- 
mation), and S will therefore contain a contribution from 
both differential Galactic rotation and the star's peculiar 
motion. Neglecting the peculiar velocity contribution for 
the moment, we can write 



S = (vj,— - ulsr ) sin^, 



(9) 



where is the average rotational velocity of the stel- 
lar population, and wlsr = 220 kms -1 . Because the OB 
stars do not follow Stromberg's asymmetric drift equation 
( Dehnen fc Binncy 1998] ), we simply take = t>LSR — 
3.0(r — r Q ), found to be c onsistent for this stellar popula- 
tion ( Drimmel et al. 1997 ). 

Finally, adopting (U,V,W) Q = (9, 5, 7) kms" 1 



(Drimmel et al. 1997), with positive U being radially in- 
wards, we can calculate an observed W for each star 
using Eq. (||). (Our adopted value for the Solar Mo- 
tion is comparable with other recent determinations (i.e. 
Dehnen & Binney 1998), but it should be noted from Eq. 
(|8|) that Wq enters only as an offset term.) The stars are 
binned in galactocentric radius to find the observed sys- 
tematic vertical velocity (W) as a function of r; this is 
shown in Fig.|ll|. The agreement of the complete portion 
of our sample with the entire sample is consistent with 
the assumption that stars fainter than 7.5 form a kine- 
matically unbiased sample. 

By ignoring the peculiar stellar velocities parallel to 
the Galactic plane, in Eq. (||), we have effectively ig- 
nored a (S 1 * tan b) contribution in each bin. This is done 
from necessity; we cannot calculate S* for stars from the 
Hipparcos catalogue as we do not have their complete 
space motion. However, from Monte Carlo tests we find 
that (S* tan 6) is typically less than 0.05 kms -1 . 

We now compare the observed vertical systematic ve- 
locities (W) with those expected from a long-lived warp. 
As the vertical velocity W of each star is with respect 
to the LSR, W = W* + v z (r,(j>) — v z (r Q , 0) , W* being the 
star's peculiar vertical motion, and V z (r, 4>) the systematic 
vertical velocity at (r, <p) resulting from the warp. Upon 
averaging (W) = v z (r, <f)) — v z (r & ,0), the difference in the 
systematic vertical motion at the star's and the Sun's po- 
sition. 

To find the systematic vertical velocity v z (r, (f>) due to 
a warp as described by Eq. (Q), we consider Jeans continu- 
ity equation of stellar dynamics. For the case of no radial 
motions we have 



dv_ d(vv z ) 1 d{w<p) _ 
dt dz r dd> 



(10) 



If the stellar density v is purely a function of z' — z — Z w 
and dv z /dz = 0, the systematic vertical motion can be 
shown to be 



The systematic motion induced by a warp is simply the 
result of the stars rotating with respect to the warp struc- 
ture. In our synthetic catalogues this systematic velocity 
is added to stars if they are in the warp (i.e. r > r w ). 

Fig. [ll] shows the systematic (W) for a typical syn- 
thetic catalogue generated with a non-precessing warp 
having the parameters obtained from the sloped stel- 
lar distribution. This demonstrates that a non-precessing 
warp, consistent with the sloped distribution of our sam- 
ple, should be easily detectable. The (W) of the synthetic 
catalogue lies systematically below the theoretical curve 
primarily due to a bias introduced by the error in the 
photometric distances. 



7. Effect of bias, amplitude, and precession 

The bias introduced by distance errors can be understood 
after some consideration. Stars whose distance is overesti- 
mated are likely to have larger distance errors than those 
whose distance are underestimated, as the distance error 
is normal in the distance modulus. Together with the fact 
that in a magnitude limited catalogue there are more stars 
closer to the Sun than further away, the skew in the dis- 
tance errors means that stars with larger observed dis- 
tances, d Ql are more likely to be stars whose true distance, 
d, is smaller. The effect of the misplacement of a star to 
larger distances is that its measured relative vertical mo- 
tion Wo (d ) will be smaller than its true relative vertical 
motion, W(d). This can be seen from Eq. (||), for if differ- 
ential Galactic rotation and vertical peculiar motions are 
neglected, we can write W(d) — (4.74^b/ cosb)d + Wq, 
and similarly for W (d ). Eliminating the common fac- 
tor (4.74/Xf,/ cos&) it can be shown that W (d ) — W(d) = 
(W(d) - W Q ){d /d - 1). Hence, as long as W(d) < Wq 
and d /d > 1, W (d ) will be less than W(d). As we 
look out from the Galactic center r is proportional to d, 
so if d > d due to error, r Q > r, and from the above 
W (r ) < W(r). Finally, the kinematic signature of a warp 
increases outward, that is W(r) < W(r Q ), allowing us to 
write H'o(ro) < W(r ); the average relative systematic 
velocity measured at a given r a is smaller than it's true 
systematic relative velocity at r Q . 

In addition to the above bias there are two other possi- 
ble reasons why the observed systematic motions may be 
lower than that predicted for a long-lived warp. The first, 
trivially, is that the actual amplitude of the warp may be 
smaller than we have estimated from the observed slope 
of the spatial distribution. This may indeed be the case as 
radio data yield a warp ampl itude between .3 and ,4kp c 
at a galactic radius of lOkpc ( Burton fc Hartmann 1994 ), 



v z (r, 4>) 



+ uj p I h(r) cos^ 



(11) 



while our warp parameters produce an amplitude of .8 kpc. 
In Fig. [l^ the signature expected from our standard warp 
(rh = 15 kpc) is shown together with the signature of a 
warp with half amplitude (rh = 30 kpc) . 

Lastly, the warp may be precessing in the direction of 
Galactic rotation. Fig. fL2 shows the expected signal of a 
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Fig. 11. The systematic vertical velocity (W), relative to the LSR, with respect to galactocentric radius. Filled circles 
are for 4250 distant Hipparcos OB stars to 13th magnitude, while open circles show the signal from the complete 
portion of our sample (894 stars to magnitude 7.5). The thick solid line is v z (r) —v z (r®) (using Eq. (pT|)), the velocity 
along the line of nodes for a long-lived non-precessing warp with an amplitude consistent with the observed spatial 
distribution, while the diamonds show the observable signature from a single simulated catalogue with the same warp 
parameters. The symbols show the median W in each radial bin, and their error bars are the standard deviation of 
the mean of W . The size of the bins increase exponentially from r© = 8 kpc with a incremental factor of e' 2 , with the 
exception of the first and last bins that are e 6 larger. 

warp of half amplitude precessing at — 13kms _1 kpc -1 . 
While these three effects, taken together, have signifi- 
cantly reduced the expected velocity signature, we still 
have not recovered the negative velocities seen in our Hip- 
parcos sample. Could higher precession rates, or higher 
errors in the photometric distances produce the velocity 
signature that we detect? The fourth model in Fig. p"2| 
shows the effect of a precession rate of —25 kins" 1 kpc - 
on a warp of half amplitude with elevated photometric 
error [cfm = 1. + max(0, (m — 6.5)/6) — M/12). In this ex- 
treme case negative vertical motions are finally achieved. 

8. Discussion 

We have confirmed with a larger and more complete 
sample of distant OB stars that the observed kinemat- 
ics and structure perpendicular to the Galactic plane are 
together inconsistent with the hypothesis of a long-lived 
non-precessing warp. However, with a revised model for 
the errors, negative systematic velocities can be produced, 
but only with excessively large precession rates, a smaller 
warp and large photometric errors. A precession rate even 



as high as — 20kms -1 kpc - is questionable on physical 
grounds, for it would require that the warp was precess- 
ing faster than Galactic rotation at a galactocentric radius 
of only 11 kpc! 

Could the photometric error parameter even be as high 
as one magnitude? The TTp-iTt distribution of models gen- 
erated with such errors (see Fig. Q), compared to the ob- 
served 7Tp-7r t distribution (Fig. ^) would suggest not. 

If we do not accept such large photometric errors and 
precession rates we are left with the same conclusions that 
we arrived at in our previous work; either the Galactic 
warp is not long-lived, or there are other systematic mo- 
tions present in the Galactic disk. Galactic companions 
could be responsible for a short-lived warp, or at least 
modify the velocities, but the regularity of the Galac- 
tic warp out to at least 16 kpc from the Galactic center 
makes this seem unlikely. Nevertheless, numerical simula- 
tions will have to be carried out to explore this possibility. 

Alternatively, the presence of other systematic vertical 
oscillations in the Galactic disk could mask the kinematic 
signature of the warp, but should also be evidenced by ac- 
companying deviations from the Galactic plane. The warp 
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Fig. 12. The expected systematic vertical velocity from simulations for four cases: a standard warp (*), a warp with 
half amplitude (x), a warp with half amplitude and a precession of — 13 kms -1 kpc -1 (o), and a warp with half 
amplitude and a precession of — 25kms _1 kpc _1 (A). In the first 3 cases the cm is as given by Eq. ([!]), while in the 
fourth model an additional 0.5 magnitude error is added. Thirty simulated catalogues were generated for each case, 
and for each catalogue median velocities were determined for each radial bin. The systematic velocity of each bin 
shown above represents the mean of the 30 median velocities, while the error bars are the standard deviation of the 
30 median velocities. The filled circles for the data and the solid line are the same as in Fig. 0, and are shown for 
comparison. 



amplitudes as inferred from the OB stars found here (.4 to 
.8 kpc at r — 10 kpc) are larger than those found in the HII 
gas (.3 to .4 kpc at r = 10 kpc, Burton & Hartmann 1994) 
and the dust (.3 kpc at r = 10 kpc, Freudcnrcich 1998f ), 
and may suggest the presence of such structures. In ad- 
dition Dehnen (1998) also finds that the kinematics of 
nearby stars are not consistent with a Galactic warp that 
begins within the Solar Circle and could thus produce a 
local tilt in the stellar distribution. Indeed, the observed 
slopes and anomalous kinematic signature may not be due 
to a warp at all, but to some other as yet unidentified phe- 
nomenon. 



The possibility of other systematic motions being 
present in the disk complicates any evaluation of the kine- 
matics of the Galactic warp. Until a deeper and more ac- 
curate astrometric survey is completed the possibility of 
other kinematic effects cannot be excluded nor, if present, 
be disentangled from the expected kinematic signature of 
the warp in the stellar distribution. Such a survey will 
eventually be given by a future astrometric mission like 
GAIA, which will also, via parallaxes, largely dispense 



with the bias introduced by the photometric distance er- 
rors. 
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